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Method and system for compensating non-linearities and time- 
varying changes of a transfer function acting on an input 

signal 



Description 

Field of the Invention 

The invention is directed to a method for compensating non- 
linearities and time-varying changes of a transfer function 
acting on an input signal. This input signal is first subjected 
to a first, adaptive transfer function and thereafter to a 
second, unknown and varying transfer function to generate the 
actual output signal. When the second transfer function varies, 
the first transfer function is updated for compensating these 
changes. The described method and system are therefore adapted 
to reduce the errors between an ideally expected output signal 
and the actually generated output signal. 



Background of the Invention 

For example, in the field of mobile communications, there is 
sometimes the need to linearise a transmitter's non-linear 
characteristic, in an adaptive manner, when the characteristic 
exhibits slow changes caused by temperature, bias, ageing or 
the like. One of the existing compensating possibilities is the 
digital predistortion of the input signal before applying same 
to the transmitter's power amplifier. Thereby, the error 
between an ideally expected output signal and the output signal 
actually generated in response to a current input signal can be 
minimised. This need for compensation not only occurs in the 
field of mobile communications but also in other fields where a 
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transfer function varying in an unknown manner is to be 
compensated by adapting an adaptive transfer- function. 

When the characteristic is changing, it takes some time until 
the system has been adapted to the new situation. It is 
desirable to perform this adaption as swiftly as possible. 
However, normally, a great number of iterations is necessary 
until the system has actually been adapted to the new 
characteristic. 



Summary of the Invention 

It is therefore an object of the present invention to provide a 
method and system which are able to swiftly adapt to a changing 
transfer function and thereby quickly minimise an error between 
the actual output signal caused by a current input signal, and 
an ideally expected output signal. 

0 The invention provides a method for compensating deviations of 
an unknown transfer function from an expected transfer 
function, and/or for compensating time-varying changes of a 
transfer function acting on an input signal, for minimising 
errors of an output signal generated in dependence on the input 
5 signal, which input signal is subjected to a first, adaptive 

transfer function and to a second, varying transfer function to 
generate the output signal, the first transfer function being 
updated for compensating deviations of the second transfer 
function from an expected transfer function and/or for 
30 compensating changes of the second transfer function, wherein, 
when updating one point of the first transfer function for a 
current input signal value, the first transfer function is also 
updated for at least one other point corresponding to a 
different input signal value. 

The invention furthermore provides a system for compensating an 
unknown and/or varying transfer function acting on an input 



35 



PCT/EP00/00579 

WO 01/56146 

signal, for minimising errors of an output signal generated in 
dependence on the input signal, which input signal is applied 
via a first, adaptive compensating means having a first, 
adaptive transfer function, to a second means having a second, 
5 unknown and/or varying transfer function to generate the actual 
output signal, the first transfer function being updated for 
correcting deviations of the second transfer function from an 
expected (ideal, wanted) transfer function by means of a 
processing means, wherein the processing means is adapted to 
10 update, when updating one point of the first transfer function 
for a current input signal value, the first transfer function 
also at at least one other point corresponding to a different 
input signal value, 

15 Due to this multi-point-correction of the first transfer 

function, the described method and system are able to quickly 
adapt the first transfer function to changes of the second 
transfer function and/or to quickly adapt the overall transfer 
function to the desired one. Thus, the apparatus swiftly 
20 converges to the new (or desired) condition, with a drastic 
reduction of the necessary adaption time. The occurrence of 
errors between an output signal ideally to be expected for a 
given input signal, and the actually generated output signal is 
therefore limited to a very short time interval after a 
25 variation (or the first use) of the second transfer function. 

Furthermore, the error deviations within this time interval are 
reduced to smaller values. 



30 3rief Description of the Drawings 



The present invention will be more readily understood upon 
referring to the following description of preferred embodiments 
when read in conjunction with the accompanying drawings. In the 
35 drawings, 

Fig. 1 shows a basic configuration of a functional system, 
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Fig. 2 shows a chart for explaining the algorithm used in the 
present invention, 

Fig. 3 shows the system of Fig. 1 in a practical implementation, 

Fig. 4 shows an iterative updating of two values of a transfer 
function, 

Fig. 5 illustrates the updating of a plurality of values of a 
transfer function, 

Fig. 6 shows an embodiment of the invention implemented as an 
adaptive predistorter and amplifier, 

Fig. 7 shows an example of the magnitude and phase response of 
function g, 

Fig. 8 illustrates an example of the ideally expected and 
actually generated output signals for a non-linearised function 



Figs. 9 and 10 show simulation results for the presented 
algorithm and known algorithms, 

Fig. 11 illustrates F values calculated at the end of the 
simulation, 

Figs. 12 and 13 show the in-phase and quadrature-phase 
components of input and output signals, 

Figs. 14 and 15 illustrate the different convergence behaviour 
of the presented and known methods for different variations of 
the transfer function g, and 

Figs. 16 and 17 show the results of three different updating 
procedures usinq the linear and LMS methods, resp. 
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Detailed Description of Preferred Embodime nts of the Invention 

Generally, the described adaptation algorithm is a method of 
minimising the error |y-x| between y=g(f(x)) and x, where g is 
an unknown function and f the adaptive function for which the 
characteristic is changing to track g. Ideally, 
f=g" J . Unlike current algorithms, it updates not only the 
transfer function f at the current input value x, but also the 
transfer function f at other points (corresponding to different 
input values x update ...) • One of the applications for such an 
algorithm is digital predistortion where a transmitter's non- 
linear characteristic needs to be linearised, in an adaptive 
manner since the characteristic exhibits slow changes with 
temperature, bias, ageing.... 

First, the new. algorithm used in the method and system will be 
mathematically described. 

Let us consider an unknown function g defined as: 
g: D g c C -» C 

x -> g(x) 

5 c designates the complex domain, and can be extended to an n- 
cimension domain; D g : definition domain of g. 

The adaptive function f is defined as: 
f : Df c C -» C D f : definition domain of f 

30 x f ( X ) C: complex numbers 

The functions (transfer functions) f and g are performed in the 
order shown in Fig. 1. Fig. 1 illustrates a functional system 
having a compensating means 1 to which an input signal x is 
35 suoplied at a first input terminal. The compensating means 1 
generates an output signal z in accordance with its transfer 
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function f (correspondence between x and z),and the value of x. 
The signal z is input into a second means 2 which outputs the 
output signal y and has a transfer function g which may change 
over time. The output signal y is fed back to a second. input 
terminal of the compensating means 1. As the overall transfer 
function is to be maintained constant, any change of the 
transfer function g or deviation from the ideal transfer 
function will be compensated by appropriately adapting the 
transfer function f so as to maintain the desired relationship 
between x and y. The means 1 may for example be an adaptive 
predistorter whereas means 2 may be a power amplifier. When the 
desired relationship between x and y is equity (amplification 
gain factor = 1) for example, any deviation of the gain of the 
amplifier 2 from the value 1 will be compensated by setting the 
amplification characteristic of means 1 to the reciprocal 
value . 

The algorithm attempts, in the present example with unity gain 
(gain factor =1), to minimise the error I y-x I between the 
20 output y and the input x, in an iterative process (see Fig. 2). 

Fig. 2 is a diagram showing the characteristic (transfer 
function) g and the desired linearised overall response 
characteristic 3 (x=y) . The horizontal axis represents the 
25 absolute value of the output z of means 1 whereas the vertical 
axis represents the final output y (absolute value) . 

m the following, f*(x) denotes f(x) at iteration k, where an 
iteration is the process of updating simultaneously one or more 
30 values defining the function f. It will also be written as: 
f k (x)=x.F k , x (with F k , x eO . 

At iteration 0 (k-0), the transfer function f is set to 1: 
VxeC, F 0 ,x=l+j-0 
35 (i.e., VxeC, f 0 (x)=x) . 
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y = g( Z ) = g(f k (x)) = g(X.F k ,x) = X.Fk.x.Gk.x 
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(1) 



_ y g(x-F k ,J 
with x.F klX =fk(x) and G Kx =--—^— 



Ideally, we would like: 
5 V xeC, y = x 



substituting in (1) , 

(2) 

o x.F k;Z .G k ,x - X 



10 oF opt =-^- (defined as such) 



(3) 



Fopt would be the optimum value to use for F k(X to satisfy (2) . 

Unlike previously published algorithms (linear method, secant 

method, Least Mean Square (LMS) etc., which have been used more 
15 specifically for predistortion systems), it is not 

theoretically correct to update the new value F WM with F opt 
(even though adaptation factors are used to avoid 
instability). If F k+1 ,x=F opt would be given, (2) would become: 
x.F k+1 ,x.G K(X = x which is the desired result. 



20 



However, this does not take into account the fact that changing 
Fw .k would mean x.F k+1 .** x.F k „, yielding a different input to 
g (see Fig. 2) . 

25 Let's define yi and y 2 as: 

y. = y =g(f k (x)) = g(x.F klX ) - x.F k , x .G k ,x 

, „,„ r \ = x Ft- v.G k+ i x (with F k +i,x=F op t here) 

y 2 = g(r k +x(x)) = g(x.F k +i,x) x. c kT .,x-v-^i,x 

therefore y 2 * x because Gfc,x * G^;, x 



30 



What really ought to be updated is the value of f corresponding 
to a different input x updale , while remains unchanged 

(F k+1 ,x=F k ,x) • x updace is defined as: 
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f" K+ l ( Xupdate) = AM . .. 

<» Xupdate- F k+ i, xupdate = X.F K ,x With F k+ i, xupdate = F 0 pt 

X-Fm (4) 



5 O Xupdate 



p 

1 opl 



substituting for F opt from (3) yields, 

Xupdate = X.F k , x -G k , x 

10 

comparing (5) with (1)/ 

O Xupdate = y 



(5) 



(6) 



Therefore, now (iteration k+1) , x up date yields the same input to 
15 g as was x at the previous iteration (iteration k) : 

/k+l( X update ) = A 00 

g(A + .(x„ pdale )) = g(/ l! W) 

■ — — v / 

assuming j; lias not varied from 
iteration k to k+l 

<^ Gk+i, xupdate = G x ,x 

O y U pdate= g(fk + l (Xupdate) ) = 9 ( ^update • F k+ 1, xupdate ) = g ( Xupdate • F 0 pt ) 

20 o y uP date= g(x.F k , x ) from (4) 

(7) 

O yupdate ~ y 

Recalling (6) and -having (7), the result is 

(8) 

25 Yupdate ~~ ^update 

which corresponds to a linearised system at point x u?da te 
(remember the aim was to have y=x) . 

30 In this form however, the algorithm may present the danger of 
not updating f completely since the current f k . x is pointing at 
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another value f MlMt »» to be updated. And once the latter has 
been updated, it will not point to another f k+x .x3 (x 3 *x up date, 
determined in the same manner as x upda t e was obtained from x) but 
merely to itself. 

Bearing this in mind, the problem is preferably overcome by 
updating the current value f M , x for the current input x too. 
The following overall solution is suggested: 



(a) 



F k+ i. xupdate = /i-Fopt or similar 
F k+ i,x = chosen algorithm (linear, LMS, secant etc..) (b) 
where //eR (real numbers) and is constant (, .adaptation factor") . 



15 



20 



25 



It is important to note that the key point here is to update 
not only the F k(X value corresponding to the current input x 
(i.e. Fj+i,x _ this is done in (b) ) , but also F k , xup date 
corresponding to a different value x updat e defined by (4) (i.e. 
F k n,xupdate - this is done in (a)). Hence, a different algorithm 
could be used to update F k(X (F k+1 ,x) , like LMS, the linear 
method, Recursive Least Squares (RLS) etc... In other words, the 
important part is (a), whereas any algorithm can be used for 
(b); the adjunction of (a) to (b) improves (b) . 

One could choose for (b) : F k+1(X = with (a,/1)eB?. 

Note that, for example, if a=/?=/<=l and assuming the algorithm 
converges, the following holds true: 



30 F^F,, 

u,a,p can be chosen as a function of the input x, and therefore 
be different for different x inputs. 
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in the foregoing explanation, the equation x=y was assumed as 
desired relationship between x and y. In practise, however, the 
present algorithm can also be applied, with appropriate 
changes, to a more complicated dependency of y on x, for 
5 instance: y-X.x (XeC, the previous description used W) . 

It can furthermore be envisaged to update adjacent indices to 
F ktl(Xupdat e as well (corresponding to other inputs) , with 
appropriate weighting.' This is because g can be assumed to have 
10 a continuous derivative (g is „ smooth-) ; doing this results in 
an even faster convergence (a particular illustration is 
discussed in the next section below) . 

One possible variation consists in updating f for all values 

u., interoolation. A possible way of doing 
15 between x and x U p da te, by interpolation. " 

this is to use linear interpolation. Assuming x < x upda te, 
P -F 

F, , • , , =F k+ls + — interpose ' 

*• update n 

With X taeniae 6 [^update J 

This can increase convergence speed quite considerably. 

in the following, a practical implementation will be described 
where D f is finite. 

in practise, the above described algorithm may be implemented 
25 in a comouter program or in a DSP (digital signal processor), 
requiring digitisation. It is then suitable to limit the srze 
of D t and define f as a list of complex values in a table „hrch 
a,e constantly updated. The table can be multi-dimensional ana 
the addressing customised to the application (e.g.: one coulo 

1-hP inout value to address the table, or 
30 use the magnitude of the input 

the square magnitude/power etc..) . 

The following example is an illustration of this particular 
case with a one-dimensional table of sire m, addressed by the 



20 



WO 01/56146 PCT/EP00/00579 

11 

magnitude of the input x. Since the table is finite, each input 
value x is associated with an index to the table 
(quantisation), fetching the corresponding value F i(k in the 
table. For example: f k (x) = x.F i(k ; f k (x) can also be produced 
5 by interpolation. 

Figure 3 shows a practical implementation of this system 
similar to same of Fig. 1 wherein the compensating means 1 is 
implemented as a table memory (here: one-dimensional look-up 
10 table) 4 which is addressed by the actual value of the input x. 
The values memorised in the table memory 4 are used for 
multiplication with the value of x for generating z, and are 
iteratively updated. 

15 Figure 4 illustrates the update process and shows, in the left 
part, the table contents (values F 1(k to F m , k ) at iterative step 
k, whereas the right part of Fig. 4 illustrates the table 
contents (values F 1(k+ i to F m , k+1 ) at the next iterative step k+1. 
The arrows between the left and right parts of Fig. 4 indicate 

20 the double updating not only of F i(k+1 but also of F jlk+ i/ 
corresponding to x and x U pdate? resp. 

Applying the interpolation discussed above results in the 
diagram shown in Fig. 5. This can increase convergence speed 
25 considerably, however, it is computationally more complex and 
resource consuming. Fig. 5 illustrates the structure wherein 
more than two F values are updated based on the present value 
of x. Similar to Fig. 4, Fig. 5 shows, in the left part, the 
table contents (values F 1(k to F,,-J at iterative step k, whereas 
30 the right part of Fig. 5 illustrates the table contents (values 
Fi, k+ i to F»,k+i) after the next iterative step k+1. Here, in 
addition to the updating not only of F i(k+ i but also of F 3 , k+ i all 
values of F assigned to values of x lying between x and x u?d ate 
are updated |P t . llW ... F j+1(k+1 ) , for example by interpolation. 
35 The plurality of arrows between the left and right half parts 
' of Fig. 5 indicate this multi-updating based on one actual x 
value . - - _ . . 
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The algorithm described above has to be adapted to this 
application, but basically remains very similar. We can choose 
for instance to address the F table by the magnitude of the 
5 input x (affects look-up table resolution) . 

It is assumed that the input x has a known limited magnitude 
range: I x I e [ I x I min ; ULxL The F table size is set accordingly 
to cover all the input values of x within this range. The total 
10 range of magnitude of input x is divided into subranges of 

possibly, but not necessarily egual size (it can for instance 
be a function of the input signal's statistical distribution 

. . . ) , the number of subranges corresponding to the number in of 

values of F contained in the table. Therefore, each actual 
15 value of x addresses that table value F which is attributed to 

that subrange to which the present x belongs, for instance F L . 

Mathematically: 

for |x| 6 Ui, a i+ i], x 3 index i in the F table corresponding to F ik 

associated to 

with aieR + , and ie[l-»m] 

Hence, f k (x) = x. F i(k (interpolation may be used here to 
determine f k (x) ) . 

The index i can be determined as: 



20 



i = floor 



(9) 



25 Equation (9) makes sure that the indices i range from 1 for 
|x|=0 to i=m for |x| = |x| ma x- 

Repeating the algorithm as it has been described in part 1: 
y = g(fY.U)) = g(x.F ilk ) = x.Fi.n.Gi.k.x 

y g(x-Fj, k ) 

30 with x.F i(k = f k (x) and G ik>x = ; — 

z x.r jk 



Fopt - — 



i.k 

1 
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As this practical implementation of the system has a finite D f , 
the additional index "i" is added which refers to the value F 
in the look-up table corresponding to cell "i" (see Fig. 4). 
»k" is the iteration number, whereas "x" is used as a referring 
5 index for the gain g at point x.F i(k which is the input to the 
transfer function g. 

So Xupdate = i^i an d the associated index in the F table for 
Xupdate is j, with j*i. 

10 

It must be ensured that |x update l is within the range covered by 
the F table (i.e., |x| rain ^ lx upda tel ^ l*lmax <=> je[l-» 
m]) . 

15 If |x updat el ^ IxUn then j is preferably set to j=l 

and if |x| raax < lx updat el then j is preferably set to j=m. 
Alternatively, the F table could simply not be updated at index 
j for jg [l->m] . 

20 Finally, the following settings are selected, similar to the 
above discussed case: 

P -a F U) 
F i(k+1 = chosen algorithm (linear, IMS, secant etc..) (ii) 

25 

v;here /i£R and is constant. 



If i=j, one may 



choose either (i) or (ii) for the update. 



30 This algorithm has been simulated in a system commonly referred 
to in the literature as predistortion. 

Fig. 6 shows a predistortion block diagram. The input x is 
complex and has been represented as: 
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x(t)=sin(co.t)+j .sin(co.t) => 



14 

|x(t)| = V2x|sln(6).t)| 
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co=2xnx8 .751e3 rad.s" 1 . 

Sampling time=iteration f requency=200e3 Hz. 
5 F table size is m=128 . 

(The expression "e3" means a multiplication by a factor 10 ). 

In Fig. 6 and the following figures, the following applies: 

10 i in =Re(x); q in =Im(x) [input - x] 

i pd =Re(z); q pd =lm(z) [predistorter out - f M ] 

i £e =Re(y); q £e =Im(y) l„™ out- - g(f(x))] 
(PA = power amplifier) • 

15 The adaptive predistorter amplifier shown in Fig. 6 comprises a 
multiplier 5 to which the input signal x and the output of a 
look-up table (memory) 6 is applied. The contents of look-up 
table 6 correspond to the table shown in Fig. 4 or 5. Thus, the 
output signal of look-up table 6 is the F i(lt value selected 
according to the present amplitude of input signal x. The 
multiplier 5 generates, as its output, the signal z which is 
applied to means 2, similar as in Figs. 1 and 3. Here, the 
means 2 is an amplifier having a time-varying gain 
characteristic. The output signal y of means 2 and the input 
signal x are supplied to a processing means 7 which stores and 
processes the above described adaptation algorithm. Updates for 
Vne indices i, j, etc., i.e. F values F lffcrt , F jf »i, etc... for 
the present iteration step "k+1" are calculated within 
orocessing means 7 based on the algorithm and are then supplied 
from processing means 7 to the look-up table 6 for storage 
therein, replacing the previous respective values (see Fig. 4 
or 5) . 

The complex input x is not only supplied to processing means 7 
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but also to a squaring means 8 which generates the square value 
of the absolute value of input x. The output of squaring means 
8 is connected to an input of look-up table 6. The present 
output value of squaring means 8 is used as address (index i) 
5 for selecting the F ilk value to be updated, so that same is 
selected depending on the square of the absolute value of x. 
This selected F value and the at least one further F value (see 
Fig. 4 or 5) are then updated according to the calculation 
results of processing means 7. 

10 

The characteristic of g (means 2) is shown in Fig. 7. The upper 
diagram of Fig. 7 represents the magnitude response of g 
(horizontal axis: absolute value of input signal z; vertical 
axis: output y- g(UI» whereas the lower diagram illustrates 
15 the phase response of g (horizontal axis: absolute value of 
input signal z; vertical axis: phase of g(|z|> in radians). 
These curves correspond to typical measured data for a given RF 
power amplifier (data has been normalised) . Typically, such an 
amplifier will produce essentially AM to AM and AM to PM 
20 distortion. 

For reference, Fig. 8 shows the input and output signals for a 
non-linearised function g (i.e., VxeC, ft(x)-x o V(i,k)e[l- 
m]x N , F ilt (x)=l +j .O =» y=g(f(x))=g(x)) with predistortion 
25 being switched off. The >u PP er diagram of Fig. 8 illustrates the 
time-behaviour of the real parts of input x (i in - Re(x>) and 
cutout v (i fe - *e(y)) with regard to the number of iterations 

(horizontal axis) . The lower diagram of Fig. 8 illustrates the 
time-behaviour of the imaginary parts of input x 1<U„ = Im(x) 
30 and output y (q £e - Im(y>> with regard to the number or 

nations (horizontal axis). Here, we are aiming at a gain 

factor of 1 so the input signals are similar to the ideally 

expected output signals. 

~. ! o^or- a case where the function g is unknown 
35 Let us first consider a case wne-~ 

T , c v^nqfer function is shown in Fig. 7. 
but time-invariant. Its transre. runo 
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for the presented algorithm. The upper diagram shows the 
imaginary part of input x over the number of iterations whereas 
the lower diagram shows the imaginary parts of signals z and y. 

5 it can be seen from the simulations that the presented 

algorithm presents faster convergence than the linear method or 
the LMS method. 

Next, let us consider a case wherein g increases monotonically 
10 with time. 

The new unknown transfer function g 1>2 is defined as a ramp of 
slope 1 multiplied by the previous transfer function g: 
gi . 2 (t,x)=g(txx) where g is the same function as defined in Fig. 
15 7. 

Again, Fig. 14 presents the results for this time varying 
function g (ramp) for all three methods, with the inscriptions 
of the horizontal and vertical axes similar to same of Figs. 9, 
20 10. 

Here again, the proposed algorithm reveals to converge faster 
and hence to be more robust against slow variations of the 
function g (gain) . 

in the next considered case, the transfer function g is a time 
varying function g in form of a step. This new transfer 
function gi . 3 , is defined as the previous function g multiplied 
by a step function: 

gi .3(t,x) = gdxx) for 0<t<0.01 (t-0.01 <=> iteration 2000) 
gi.3(t,x) = g(1.2xx) for 0.01<t 



25 



30 



35 



Fig. 15 presents the results for all three methods 
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This case attempts to test the response of the algorithm for 
sudden variations of g. The proposed method again converges the 
quickest. Note that the LMS method is quite slow to readapt. 

A comparison for g as a step (same as above) is done for the 
linear method and the LMS method for the three following 
algorithms (Figs. 16 and 17): 

1) Updating current x address only (index i) ; 

2) Updating current x (index i) and x upd ate (index j) 
addresses; 

3) Updating current x (index i), x upda t e (index j) and all 
xinterpoiate addresses between x and x upda te (indices w with we 
[j->i] if j<i for example) . 

The algorithm 1) is the known one whereas algorithms 2) and 3) 
are in accordance with the present invention. 



Fig. 16 shows the three updating procedures using the linear 
0 method, whereas Fig. 17 shows the three updating procedures 
using the LMS method. The upper diagram illustrates algorithm 
1), the middle diagram algorithm 2), and the lower diagram 
algorithm 3) . 



>5 



30 



It is obvious that algorithm 2) has a better convergence speed 
than algorithm 1), and that algorithm 3) yet improves this 
speed dramatically. Interpolation used in algorithm 3) relies 
on the assumption of a „well-behaved" response for g, which 
becomes even more plausible as the algorithm converges. In 
general, for a power amplifier, this is a fair assumption. 



The algorithm described here thus provides, compared to other 
algorithms used for predistortion systems such as the linear 
method or the LMS method, very fast convergence, and is quite 
35 robust against fluctuations of g in time. Besides, it requires 
few operations (i.e. is computationally simple), and can be 
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used in conjunction with other algorithms as an enhancement. 

The use of spectrally more efficient modulation schemes (higher 
throughput/occupied bandwidth ratio) to accommodate higher data 
rates in communications often introduces both phase and 
.agnitude variations of the signal, thus requiring stringent 
linearity of the transceiver to preserve signal integrity. 
Transmitter linearisation, in particular, is then often 
necessary to provide good linearity while offering reasonable 
power efficiency. 

Satellite and mobile applications are two potential candidates 
for linearisation. The proposed algorithm can, for instance, be 
used within digital predistortion systems as a linearisatxon 
means. The present invention can therefore be embodied in a 
stationary or mobile device £or mobile oommunications xn 
particular in a RF transmitter thereof, in such standards as 
GSM EDGE (Global System for Mobile Communications Enhanced Data 
Rates for GSM Evolution) , TETRA (Terrestrial Trun(ced Radxo, 
CDMA (Code Division Multiple Access,, W-CDMA (Mi deband-CDMA> 
IMC Worth American Digital Cellular), PDC (Personal Dxgxtal 
Cellular) . . . 

An embodiment is therefore a stationary or mobile 
, m0 bile communications containing a device as shown xn « 3 or 
6 , orefexably as defined in the system claim, and/or operatxng 
as defined in the independent method claxms. 
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Claims 



1. Method for compensating deviations of an unknown transfer 
function from an expected transfer function, and/or for 
compensating time-varying changes of a transfer function acting 
on an input signal, for minimising errors of an output signal 
generated in dependence on the input signal, which input signal 
is subjected to a first, adaptive transfer function and to a 
second, unknown and/or varying transfer function to generate 
the output signal, the first transfer function being updated 
for compensating deviations or changes of the second transfer 
function, 

wherein, when updating one point of the first transfer function 
for a current input signal value, the first transfer function 
is also updated for at least one other point corresponding to a 
different input signal value. 



3 2. Method according to claim 1, wherein the at least one other 
point of the first transfer function is calculated in 
dependence on the current input signal value and the value of 
the first transfer function at this input signal value. 



3. Method according to claim 1 or 2, wherein the at least one 
other point of the first transfer function is calculated in 
additional dependence on the respective value of the second 
transfer function. 



30 



4 Method according to claim 1, 2, or 3, wherein the at least 
one other point of the first transfer function is corrected in 
dependence on the value of the second transfer function for the 
35 current input signal value. 
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5. Method according to claim 4, wherein the at least one other 
point of the first transfer function is corrected by 
multiplying the value of the second transfer function for the 
5 current input signal value, with a factor. 



10 



6. Method according to any one of the preceding claims, wherein 
values of the first transfer function lying adjacent to the one 
point of the first transfer function and/or the at least one 
other point of the first transfer function are also corrected 
with appropriate weighting. 



15 7. Method according to claim 6, wherein values of the first 
transfer function lying between the one and the at least one 
other point of the first transfer function are corrected by 
interpolating . 
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8. Method according to any one of the preceding claims, wherein 
the values of the first transfer function are provided in form 
of a table being addressed by the input signal values. 



9. Method according to any one of the preceding claims, wherein 
the updating of the first transfer functio n is done in an 
iterative manner. 



10. System for compensating deviations of an unknown transfer 
function from an expected transfer function, and/or for 
compensating time-varying changes of a transfer function acting 
on an input signal, for minimising errors of an output signal 
35 generated in dependence on the input signal, which input signal 
is applied via a first, adaptive compensating means having a 
first, adaptive transfer function, to a second means having a 
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second, unknown and/or varying transfer function to generate 

. , ,„„ first transfer function being 
the actual output signal, the first trans 

0 pdated for compensating deviations and/or changes of 

second transfer function by means of = processing means, 

wherein the processing means is adapted to update, when 

gating one point of the first transfer 

input signal value, the first transfer function also at at 

least one other point corresponding to a different input signai 

value. 



U. use of the method according to anyone of claims 1 to 9 or 
o£ the system according to ciaim 10, to perform linearisation 
of a signal path. 
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Figure 2: Illustration of the algorithm 
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Figure 4: F table 
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Figure 5: Variation for F table update 
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Figure 6: Predistortion block diagram 
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Magnitude response of g 




Figure 7: Magnitude and phase response of g 
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Non-linearised function g- (predistortion off) 




number of iterations 
Figure 8: Non-linearised function g - (predistortion off) 



H=1; a=0.1; p=0.9 




Figure 9: 20*logl0(|y-x|) for the presented algorithm 
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Linear method - adaptation factor=0.9 




LMS method - adaptation factor=Q.9 




Figure 10: 20*logl0(|y-x|) for linear and LMS methods 
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Figure 
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13: Quadrature components tnrougn 
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Figure 14: time varying function g (ramp) 
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Figure 15: time varying 
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